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ABSTRACT 

This contribution to the Proceedings is based on the talk given at the Conference on 
Birth of the Universe and Fundamental Physics, Rome, May 18-21, 1994. Some selected 
topics of the subject are reviewed: Models of Primordial Fluctuations; Reconstruction 
of the Cosmological Density Probability Distribution Function (PDF) from Cumulants; 
PDFs in the Zel'dovich Approximation and from Summarizing Perturbation Series; 
Fitting by the Log-normal Distribution. 
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1 Introduction 



The canonical model for the formation of large-scale structure is based on gravitational instability 
of small initial fluctuations in the matter density 6. The primordial fluctuations are assumed to be 
a random field, which is fully specified by the joint probability density functions (hereafter PDFs). 
These are assumed to originate from quantum fluctuations of the inflaton scalar field Scpk^^^^ that 
were stretched to large comoving scales during the inflation phase For the simplest case of a 
single inflaton scalar field, curvature (adiabatic) density perturbations are generated by the inflaton 
fluctuations, 6 oc 50, which is a superposition of spatial waves e*'^^ with random phases, realizes a 
random Gaussian field. 

Theoretically, there are scenarios based on non-Gaussian initial fluctuations. Althogh these 
models currently do not appear to be our first choice, they are interesting as alternatives. I will 
discuss these models in the next Section. The statistical nature of the initial fluctuations is therefore 
a basic distinguishing feature between competing theories. 

The COBE DMR temperature anysotropy fluctuations sky map, in principle, can directly probe 
the statistics of the primordial fluctuations on large scales. It is consistent with Gaussian fluctua- 
tions, but its signal to noise level so far is sufficient only to rule out strongly non Gaussian models 

i- 

In any case, at present the smoothed galaxy distribution is not a Gaussian random field, due 
to the nonlinear dynamics of gravitational instability. During linear evolution, when all Fourier 
components 6k evolve at the same rate, the cosmic PDFs do not change form. However, even 
quasilinear evolution, which makes phases correlated, introduces strong non-Gaussian features. 

The study of the PDFs of the large scale cosmic density and velocity fields, and their moments 
drew much attention recently. In Section 3, I will discuss recent analytic results on the derivation of 
the density and velocity PDFs from gravitational dynamics. There were suggestions to design the 
density PDF phenomenologically ^ . For instance, the log- normal form is a surprisingly good fit 
for the densuty PDF in CDM model. Why it may be so, I will discuss in Section 4. 

The discrete analogy of the one point PDF is the counts in cell statistics. The density PDF is 
obtained by the smoothing of the discrete galaxy field with a sharp filter. The observed density 
PDF manifests significant non-Gaussian features in the non-linear and even in the quasi-linear 
regimes, and can be well fitted by the log-normal distribution. Assuming that galaxies trace mass, 
observed galaxy distribution is consisitent with Gaussian initial fluctuations on scales of these 
surveys. However, the discriminatory power of the observed PDFs is currently limited for several 
reasons. I will discuss these issues in the Conclusion. 



1 



2 Models of Primordial Fluctuations 



It has been shown that in the framework of the the inflation picture there is still room for non- 
Gaussian fluctuations. In case of the curvature perturbations it requires specific forms of the inflaton 
potential or coupling with other fields, plus fine tuning of the parameters 0. On the contrary, in 
case of the isocurvature perturbations generation of non Gaussian fluctuations is rather typical 
[^, 0] . Let another light scalar field x be present at the inflation. Long wavelengths fluctuations 5x 
inevitably arise during inflation. If later the x-particles become dominant and responsible for dark 
matter, then we come to the model with isocurvature perturbations. The fluctuations 5x are Gaus- 
sian. However, they can generate the isocurvature energy density perturbations 5isocurv = F{5x), 
which are non Gaussian for nonlinear functions F (local non Gaussian field P]). An interesting 
example is the cosmic axion as x-field. The axion energy density is proportional to the axion 
potential V{x) = ~ cosx/xo); therefore isocurvature fluctuations in the model with axions 
as CDM are non Gaussian. The power spectrum of isocurvature fluctuations does not depend 
on their statistics. Unfortunately, the observed power spectrum ranging from the horizon to the 
galaxy clustering scales leaves small room for the CDM scenario with mixture of curvature and 
isocurvature fluctuations. 



Another possibility is that the x-particles are underdominant but their fluctuations are signif- 
icant and play the role of the seeds for the structure formation . There are models of the local 
non Gaussian isocurvature baryon fluctuations §]. Non-Gaussian density fluctuations also arise in 
scenarios, where they originate from topological defects, such as cosmic strings or textures |jl3[, 
or late phase transition [|l^] or from non-gravitational cosmic explosions [0. Initial shape P{S) 
in scenarios with topological defects contains the long tail in the overdense region, for instance, 
exponential for textures. It is interesting that the velocity PDF in these model is normal, due to 
the central limit theorem 1113 



3 Reconstruction of PDFs from Cumulants 



First, let us recall basic terms which are adopted in the literature. We use the cumulants of the 
distribution < 6^ >c rather than its moments < 6^ >. The generating function of the cumulants 
is C(/i) = I]^2 < >c f^^/p^-, where /z is an auxiliary parameter. All cumulants of the normal 
distribution are vanishing, except the linear variance of the density fluctuations a. It is convenient 
to use rescaled cumulants Sp{a) =< >c /a'^^P~^\ as a set of descriptive constants of a distribution. 
The Sp parameters would be constants in the hierarchical ansatz, but in general they are functions 
of a. The parameter S3 multiplied by a is the skewness, and 6*4 multiplied by is the kurtosis. 
The advantage to use these parameters is that they are final numbers Sp{0) for an arbitrarily small 
a. The parameters Sp{0) are the fingerprints of the particular dynamical model. 

When the whole series of the cumulants is known it is possible to reconstruct the density PDF 
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from the generating function of the cumulants 



1 f+ioo 

P{p) = — rf/iexp[-C(/i)-(p-l)/x]. (1) 

ZTTZ J-ioo 

Here and further p is the density in units of mean density. In practice we may have only a few 
lower cumulants, such as the skewness and the kurtosis. Theoretically, they can be calculated in 
the perturbation theory. In the case of the weakly non-linear dynamics, when a slight departure 
from the initial Gaussian distribution is expected, one can obtain the general decomposition series 
around the Gaussian PDF induced by the first non-zero cumulants - the Edgeworth expansion. In 



cosmological context this was implicated in [14, 15, 16 . 



The Edgeworth expansion for the cosmic density PDF can be obtained from the reconstruction 
formula (|1|), keeping a few lowest terms in the generating function The result is 



Pi5) 



(27ra2)i/2 



2a^ 



Q / Q C2 



1 + a-^Hs {v) + a' [^-^^H, {v) + -^H, {v) ) + ... 



(2) 



where -ff„(z/) are the Hermite polynomials, v = bjo . The actual forms of the parameters = Sp{cr) 
which depend on the particular dynamics, affect the expansion with respect to a. The usual 
measurements of the lowest cumulants are significantly affected by the high density tail of the 
PDF, i.e. the rare events. It is interesting, that the lowest cumulants alone are responsible for the 
shift of the peak of P{S). Therefore the measurement of the shape of the PDF maximum, which 
statistically is more robust, can provide an alternative method of evaluation of the lowest cumulants. 
The Edgeworth expansion fails to reproduce P{S) when \6\ > 0.5 where spurious wiggles appear. In 
practice it is useless for a > 0.3. The Edgeworth expansion also can be used for modelling slightly 
non Gaussian initial fluctuations. 



4 PDFs in the Zel'dovich Approximation 

Physical meaning of the density PDf is the fraction of volume with a given level of density. The 
former shape of the initial density PDF is broken as non-linearity develops, since matter is evac- 
uated from the underdense regions, which occupy the largest fraction of volume, meanwhile the 
overdence regions collapse. In the formal expression (j^) this corresponds to the positive sign of the 
skewness, increasing with time. This feature of the density field evolution is clearly manifested in 
the Zel'dovich approximation 



|(i-moi)(i-mo2)(i-mo3)|' 

where po is the initial local density, growing mode of adiabatic perturbations is D{t), and Aoi are 
the eigenvalues of the Lagrangian deformation tensor. This formula does not assume any particular 
initial statistics, whether Gaussian or not. The gravitational clustering at sufficiently large scales 
can be considered in the quasilinear theory in a single stream regime ignoring small scale details. 
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The Zel'dovich approximation (j^) can be applied for the filtered initial gravitational potential. This 
approach sometimes is reffered to as the truncated Zel'dovich approximation. 

In the truncated Zel'dovich approximation the statistics of the evolved density field can then 
entirely be obtained from the statistics of the initial local density po and the initial eigenvalues Aq-s. 
For adiabatical perturbations po is the background mean density. Let the initial joint PDF of all 
involved cosmic fields, including velocity and gravitational potential, be Wo{po, Aqi, A02, A03, Uq, $o)- 



The density PDF can be obtained in the general case of arbitrary initial statistics simply as an 
integral over all involved variables except density []T3|, |TB 



Pip) = J ^[Pl(l - ^Aoi)(l - DXo2)il - DXos)\ - l]WodXoidXo2dXo3d\od^, 



(4) 



For the Gaussian initial conditions we can omit Uq and $0 iii the initial joint PDF and write it 
as Wo(po5 Aoi, A02, A03) = S{po — l)Mo(Aoi, A02, A03). The second factor is the well known joint 
distribution function of the eigenvalues of the initial deformation tensor for an initial Gaussian 
displacement field fl^. Substituting this form of Wq into the integral (H), after some tedious 
algebra, we can reduce it to the one- dimensional integral |TB|, [13 : 



Pip) 



N 
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3(£)i/3 
^ p' 



ds e 



-(s-3)V2o-^ 



+ e 



(5) 



9-53/2 



In 



and prefactor is — 
H). The density PDF calculated 



where = s ■ 5^/^ (| + cos |(n - l)7r + | arccos - 1 

the limit of small a the expression (H) transfers into the form 
numerically from formula (^) is plotted in Fig. 1. Without final smoothing, this PDF does not 
depend on the power spectrum n. The quality of approximation of actual PDF by the formula 
is increasing as n is decreasing n —>■ —3. The ra-dependence can be taken into account in the limit 
of small cr . 



For the slightly non Gaussian initial fluctuations, one can expand the joint distribution Wo{po, Aqi, A02, 
around its Gaussian form, discussed above. For this purpose it is convenient to return to six com- 
ponents of the deformation tensor, and apply the generalized Edgeworth expansion for several 
variables. From (^) one can obtain the density PDF evolving in time from slightly non Gaussian 
fluctuations in form of series around distribution (|^). One of the lesson from this exercise is that, 
in principle, the final density statistics depends on the statistics of all fields involved in the joint 
distribution Woipo, Aqi, A02, Aqs, uo,^o)- 



^03) 



5 PDFs from Summarizing Perturbation Series 

Using the perturbation theory for Gaussian initial conditions, it is possible to derive the cumulants 
of the density PDF in the small a limit. The basic assumption here is that the gravitational 
clustering at sufficiently large scales can be considered in the single stream regime. Derivation 
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of the lowest cumulants in the quasi-hnear dynamics in the single stream regime was extensively 
elaborated [|^, ^ |2^, 1^. Bernardeau ^ found an elegant method to derive the closed 



form for the generating function of the cumulants C(/i) in the limit of small o", which allows to 
obtain all parameters >S'p(0) and reconstruct the density PDF in this regime. 

The chain of equations for density perturbation S^^^ in each order p cannot be resolved. However, 
one can define the connected averages: Up =< S^^\S^^^y >c /cr'^^, and construct their generating 
function G{t) = Y1,%=i{^'tY^p/p'^-i where r is an auxiliary variable. It is remarkable that there is a 
single equation for G{t) which corresponds to the density contrast of the spherical collapsing linear 
overdensity r p^. An approximated analytical solution for this function is 



GO 



1 



- 1. 



(1 + r/a) 

For the three dimensional single stream cosmological gravitational instability a ^ 1.5 p6[ . 



(6) 



This method and results are rather general in the limit of small a. For instance [0, for the 
Zel'dovich dynamics in the one dimensional case, where the Zel'dovich ansatz is an exact solution, 
the formula (^) is valid with a = 1. For the Zel'dovich approximation in the two and three 
dimensional cases, a = N , where is the space dimension. In the formal limit ^ 1 we get 
G{t) = exp(— r) — 1, which coincides with that of the log-normal distribution. 



The transition from G{t) to the generating function of the cumulants G{fi) is given by the 
Legendre transform with the variable of the transform equal to unity afterwards ||32|, Substituing 
the solution (|^) into (|l|), one reconstructs the density PDF. Without final smoothing, this PDF 
depends on a only and not on the power spectrum index n. However, for a practical purpose it is 
necessary to filter the evolved density field. After filtering, density cumulants and PDF depend on 
the shape of the power spectrum [^, Q . Bernardeau demonstrated |]2^ , that the effect of filtering 
reduces to the a simple transformation of the generating function G(r) — > G^{t). The resulting 
shape of the density PDF, based on the filtered generating function G^ {t), is presented in Fig. 1, 
2. 



6 Fitting By The Log-Normal Distribution 



As it was noted a long time ago by Hubble the galaxy count distribution in angular cells 
on the sky might be well described by a log-normal distribution. The recent study of the spatial 
distribution of galaxy count are in good agreement with the log- normal fit [ 
normal density distribution reads 



g The log- 



P{p) 



2\l/2 



exp 



2a} 



P 



(7) 



It was found 120] that the log-normal distribution is an excellent approximation to the density 



PDF from N-body CDM simulation for the tested values of 0.3 < a < 1.5 in the tested range of 



5 



p < 10, see Fig. 1. Such a remarkable fitting inspires the thought that there might be a dynamical 
reason to manifest the log-normal features of the density PDF. In the log-normal mapping of 
the linear density field was suggested to describe its non-linear evolution. This log-normal model is 



universal for any spectral index n. Unfortunately the log-normal mapping does not work |3^. At 
Fig. 2 the log-normal distribution is compared with the analytically derived density PDF The 
log- normal shape fits an actual PDF not for arbitrary n, but for n —1. It is easy to understand, 
considering the cumulants of the log- normal distribution. For instance, we have 5*3°^ (a) = 3 + a"^ 
for an arbitrary a. 

The actual parameter is 5*3(0) = y — + 3). The parameter of the log-normal distribution 
S^^^{a) fits an actual parameter not for arbitrary n and a, but along the curve n = — 1.14 — a^. The 
cr-correction to the 6*3(0) changes slightly the numerical prefactor of a. The parameter S^^^{a) fits 
an actual parameter 5*4(0") approximately along the same curve. Thus, the log-normal distribution 
fits well in the particular region of the parameter space (n, a) around the "log-normal" curve n((T), 
but worsening outside of this region. The actual density PDF can be further approximated by the 



Edgeworth expansion around the log- normal form |35]. By chance, the popular CDM model at 
moderate a crosses this region. 

Equipped with this method, one can understand why and where the other fits such as "ther- 
modinamic" [i, EH] or the negative binomial distribution ||32|, |3^ are applicable. 



7 Discussion 



In general case, S*p-parameters as functions of the linear density variance can be presented as series 

Sp{a,n) = Sp{{),n)+Tp{n)a^ + ... (8) 



There is derivation of all coefficients S*p(0, n) p5| 



but a little is known about Tp{n) 



which describe 



a lowest (T-correction. Apparently, Tp(n) depends on n. Some N-body simulations indicate the 
tendency Sp{a,n) to increase with a faster as n decreases [Q. Parameters Sp slowly increase with 

and for CDM model 



a at quasilinear regime, Tp{n) -C 5*p(0), at least for — 1 < n < ||14 



T^. A 



possible explanation is that the spherical collapse is the dominanted form for these spectra, contrary 
to the pancakes for n —3. Another choice of the set of descriptive parameters is < 5^ > /crnt~'^\ 
where nonlinear variance is (T„/((t). It was shown in |¥T] numerically for CDM model that these 
combinations remain constants equal to Sp{0) for a wide range of (T„;. Therefore Bernardeau's 
method, formally applicable for small a, can be extrapolated across the whole quasilinear regime, 
which makes it very useful for practical application. The cx-dependence in Sp can be more significant 
for n —3. In this case the accuracy of PDF based on the Zel'dovich approximation, is increasing. 
Another related question is the hierarchical ansatz in the highly nonlinear regime. It assumes that 
Sp{a) are saturated as a 3> 1. This is still an unanswered question, however, there are interesting 



results in this regime [39, 40 . 



In practice, galaxy count in cell distributions, and consequently, the density PDFs have been 
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measured in various catalogs in many works, for example see |3T 
the IRAS surveys 



1^ for the CfA survey 

H 



33, m 



m 



for Sp parameters 



^ for SSRS survey, in 43 for clusters of galaxies, 
from the APM survey, and references therein. These data are compatible with Gaussian initial 
fluctuations, and can be fitted by the log-normal distribution, see, e.g.. Fig. 3. Thus, some "log- 
normal" features of the observed density PDF would mean that the realistic cosmological model 
has the n{a) dependence close to the CDM-type models. Unfortunately, the error bars of the 
observed density PDF increase at large p ~ 2 — 3, and significantly affected by the depth of the 



sample p! 



Therefore currently one can rule out only strongly non Gaussian models. Another 



potential problem is galaxy bias, which affects the observed deviation from Gaussianity The 



linear bias at filtering scales would preserve the reported analytic results. However, the departure 
from this simple rule is expected, as predicted by numerical simulations [H6|. 



The density PDF deviates from a normal distribution very rapidly. On the contrary, the velocity 
PDF departures from its initial distribution very slowly [[19| , The observed velosity PDF is 
shown at Fig. 3. 
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Fig. 1 The density PDF obtained by different methods for CDM initial conditions [Q. The 
points are measured in a numerical simulation at the time corresponding to the present epoch, at 
two different smoothing radii Rq = 5h~^Mpc (triangles) and Rq = 15/i^^Mpc (circles). The rms 
density fluctuation are then respectively cr = 1.52 and cr = 0.47. The error bars have been obtained 
by dividing the sample into eight subsamples. The solid line is the prediction of the Bernardeau's 
method when the smoothing effects are taken into account. The dashed line is the prediction (^ 
from the ZA and the long dashed line is the lognormal distribution (j^). 



Fig. 2 The dynamically motivated density PDF for different n, a = 0.5, versus log-normal 
distribution The solid line for n = +1,0,— 1 is the density PDF obtained by the Bernardeau 

-2 is the PDF in the Zel'dovich approxi- 



^. The solid line for n = +1, 0, 

method with the final smoothing. The solid line for n = 
mation. The dashed line is the log-normal distribution. 



Fig. 3 PDFs for IRAS 1.9Jy density and velocity fields in a sphere of radius SOh'^Mpc, after 
Gaussian filtering with Rs = 6h~^Mpc. Dashed and long dashed curves are the Gaussian and 
lognormal distributions with the same a. Also shown the errors associated with the limited volume 
sampled p!9 |. 
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